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ABSTRACT 


Optical  measurements  of  the  structure  parameters  of  temperature 
2  2 

(C,j,  ) ,  humidity  (C^  )  and  temperature-humidity  (C^q)  were  attempted  at 

0.94  ym,  10.7  ym  (CO  laser)  and  337  ym  (HCN  laser),  but  owing  to  various 

1  2 

instrumental  problems  no  measurement  of  C^  and  C^^  could  be  achieved. 

In  situ  measurements  of  the  same  parameters  were  more  successful.  Two 

techniques  of  measuring  structure  parameters  are  compared:  by  using 

spaced  sensors,  and  by  time-lagging  one  sensor.  In  this  way  the  ex- 

pected  r  dependency  of  the  structure  parameter  could  be  checked  (r  = 

distance  between  sensors,  either  by  spacing  or  by  time-lagging).  From 

the  connection  between  structure  parameters  and  spectra,  evidence  of  an 

inertial  subrange  for  the  temperature-humidity  cospectrum  is  found. 

Special  attention  is  paid  to  the  use  of  Ly-a  hygrometry  for  measuring 

structure  parameters  involving  humidity.  Finally,  relations  between 

2  2 

the  structure  parameters  C^  and  C^  ,  and  the  fluxes  of  heat  and  mois¬ 
ture  are  discussed;  an  effect  of  atmospheric  stability  was  not  detectable 
for  -z/L  >  0.02. 


1.  INTRODUCTION 

The  effect  of  atmospheric  turbulence  on  propagation  of  electromagnetic 
radiation  in  the  atmosphere  has  received  renewed  interest  in  the  last  two  decades, 
mainly  because  of  its  degrading  effect  on  optical  communications.  For  the  meteo¬ 
rologist,  however,  this  very  degradation  provides  a  means  for  studying  atmospheric 
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turbulence.  In  particular,  by  observing  the  scintillation  of  a  visible  light  beam 

2 

one  can  infer  the  temperature  structure  parameter  C  ,  defined  by 


<<TA  -  V  > 


where  -  Tg  is  the  temperature  difference  between  two  close  locations  A  and  B, 

and  r . „  is  the  distance  A-B.  The  brackets  indicate  an  ensemble  average.  The 
AB 

locations  A  and  B  should  be  so  close  that  any  temperature  difference  can  result 

only  from  inertial- subrange  turbulent  structures  (Wyngaard,  1973).  The  quantity 

<(T.  -  T„)  >  is  called  the  temperature  structure  function,  and  denoted  by  D  (r) . 
Ad  -L 

The  temperature  structure  parameter  is  one  of  the  parameters  describing  the 
state  of  the  atmospheric  boundary  layer.  Its  magnitude  is  related  to  the  sensible 
heat  flux  at  the  Earth's  surface.  In  the  case  of  a  very  unstable  atmosphere,  the 
relation  has  the  simple  form  of 

H  =  0.806  P  cpz  (Ct2)3M  (Is)172,  (2) 

-2  -3 

where  H  is  the  sensible  heat  flux  (W  m  ),  p  the  air  density  (kg  m  ),  c  the 

-1  -1  " 

specific  heat  of  air  at  constant  pressure  (J  kg  K  ),  z  the  height  above  the 

_2 

surface  (m) ,  k  the  von  Karman  constant,  g  the  acceleration  of  gravity  (ms  )  and 

2  2  -2/3 

T  the  absolute  temperature  (K) .  CT  has  dimensions  K  m  .  For  less  unstable 
situations,  the  relation  is  more  complex  and  involves  a  stability  parameter 
(Wyngaard  and  Clifford,  1978). 

Similar  to  the  above  expression,  a  relation  has  been  proposed  connecting  the 
humidity  structure  parameter  to  the  evaporation  or  latent  heat  flux: 


LEo  =  0.806  L  z  (Ct2)1/4  (Cq2)1/2  (~V/2  , 


-2  -1 
where  LE  is  the  latent  heat  flux  (W  m  ) ,  L  the  latent  heat  of  evaporation  (J  kg  ) 

2  ° 

and  Cq  the  humidity  structure  parameter  defined  as 

„  2  <(QA  -  V2> 


2 


In  this  report  we  express  the  humidity  Q  in  kg  m 

.2-6  -2/3 

kg  m  m 


-3 


so  Cr 


has  dimension 


Relation  (2)  has  been  explored  in  the  past  to  some  extent  (e.g.,  Wyngaard 
et  al. ,  1978;  Davidson  et  al.,  1978;  Coulter  and  Wesely,  1980),  but  experimental 
data  to  test  relation  (3)  are  more  scarce  (Farrall  et  al.,  1980). 


The  link  between  structure  parameters  and  optical  propagation  becomes  obvious 
if  one  considers  the  expression  for  the  scintillation  of  a  spherical  wave  of  light 
along  a  horizontal  path  (e.g.,  Wang  et  al.,  1978): 


=  0.124  k 


7/6 


L11/6 


(5) 


At2  Ct2/<T>2  +  Aq2  Cq2/<Q>2  +  2ATAQ  Ctq/<T><Q>  . 


(6) 


The  relation  (5)  is  valid  for  the  case  of  a  spherical  wavefront,  infinitesimally 

2 

small  receiver  aperture,  and  inertial  subrange  fluctuations.  0  is  the  variance 

X 

of  the  log-amplitude  of  the  received  light,  k  is  the  wavenumber  (k  =  2tt /A)  ex¬ 
pressed  in  id  \  L  the  path  length  between  transmitter  and  receiver  (m) .  The  coef¬ 
ficients  At  and  Aq  depend  on  wavelength,  humidity,  temperature,  and  pressure,  and 
C^q  is  the  temperature-humidity  structure  parameter  defined  by 


CTQ 


<(TA  ' 


v<qa  - 


2/3 

r 

AB 


(7) 


2 

is  the  refractive  index  structure  parameter.  The  first  term  on  the  right- 
hand  side  of  (6)  gives  the  effect  of  temperature  fluctuations  on  the  scintilla¬ 
tion,  the  second  term  that  of  humidity  fluctuations  and  the  third  term  the  effect 
of  correlated  temperature-humidity  fluctuations.  The  first  two  terms  are  always 
positive,  whereas  the  third  term  can  have  either  sign. 


The  quantities  A^  and  Aq  can  be  calculated  from  a  knowledge  of  air  pressure, 
temperature,  humidity,  and  the  line  strengths  and  shapes  of  water  vapor  absorption 
lines.  The  effect  of  other  gases,  like  CC^,  may  usually  be  neglected.  An  ex¬ 
tensive  calculation  of  A^,  and  Aq,  and  of  other  quantities  referring  to  optical 
propagation  has  been  performed  by  Hill  et  al.  (1980)  for  the  wavelength  interval 
4.7  ym  to  radio  wavelengths. 


**«#«***£* 
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For  the  three  wavelengths  of  interest  (0.94  yin,  10.7  yin,  and  337  ym)  it  is 
useful  to  have  analytic  expressions  for  A^,  and  so  that  their  values  can  be 
obtained  easily  for  the  given  mean  temperature,  humidity,  and  pressure  during  a 
measurement.  These  quantities  were  calculated  (by  R.J.  Hill  and  J.T.  Priestley) 
over  a  range  of  temperatures  from  -10°C  to  38°C,  for  relative  humidities  from 
0%  to  100%,  and  pressures  810  to  1060  mbar.  The  resulting  calculated  values  were 
fitted  to  simple  functional  forms  (by  R.S.  Lawrence).  The  values  of  A^,  and  A^ 
are  essentially  the  same  at  0.94  ym  as  for  visible  light;  thus 


AT(0.94  ym)  =  -77.49  PT  x  10 
and  A^(0.94  urn)  =  -59.00  Q  x  10  8  , 

where  P  is  total  atmospheric  pressure  in  mbar,  T  is  temperature  in  Kelvin,  and 

_3 

Q  is  the  water  vapor  density  in  kg  m  .  For  10.7  ym  the  functional  fit  gave 

At(10.7  ym)  =  Ax(0.94  ym)  -  2.5  x  10"23  (Tc  +  91)7’487  R 

and  Aq (10. 7  ym)  =  AQ(10.7  ym)  -  4.172  x  10~19  (T£  +  86)6'22  R  , 

where  T^  is  temperature  in  Celsius  and  R  is  the  fractional  relative  humidity.  For 
337  ym  the  functional  fit  gave 

At(337  ym)  =  AT(0.94  ym)  -  (34.54  +  2.325  Tc  +  0.0626  +  0.001527  T,3)  R  x  10~6 

and 

2  3  6 

A  (337  ym)  =  Art(0.94  ym)  +  (35.26  +  2.353  T  +  0.06247  T  +  0.00151  T  )  R  x  10 
Q  Q  K  V  c  c  c 


For  reasons  of  convenience,  we  rewrite  (5)  and  (6)  as 

at/<t>  c, 


a2  -  .124  k7/6  L11/6 

*  <Q>2 


A^/<T>2  C2 


2  2  2 

L  V<Q>  CQ 


+  1  +  2 


TQ 

V<q>  CQ  j 


(8) 


and  illustrate  with  an  example  how  temperature,  humidity,  and  their  correlation 

-2  -3 

contribute  to  the  scintillation.  We  put  <T>  =  298  K,  <Q>  =10  kg  m  (which 

2  2  -2/3 

implies  a  relative  humidity  of  40%),  C  =  0.367  K  m  (corresponding,  using 

^  9  9  _  l. 

(2)  with  z  =  4  m,  to  a  sensible  heat  flux  of  200  W  m  ),  and  C  =  0.367  x  10 

2  -6  _  9/3  _9  ^ 

kg  m  m  (a  latent  heat  flux  of  about  400  W  m  ,  using  (3)).  We  consider 

three  wavelengths:  0.94  pm,  10.7  pm,  and  337  pm.  Finally,  we  assume  a  perfect 

positive  correlation  between  temperature  and  humidity  fluctuations,  consequently 
2  2  1/2 

Ctq  =  (CT  CQ~)  .  The  values  of  the  three  quantities  between  the  square  brackets 
in  (8)  are  given  in  Table  1.  Only  for  the  longest  wavelength  considered  (337  pm), 
is  the  effect  of  temperature  fluctuations  small,  whereas  at  10.7  pm  both  tempera¬ 
ture  and  temperature-humidity  effects  have  importance  and  at  0.94  pm  temperature 
fluctuations  are  dominant.  This  observation  leads  to  the  original  goal  of  the 

present  work:  by  measuring  the  scintillation  at  these  three  wavelengths  it  should 

2  2 

be  possible  to  obtain  three  equations  with  three  unknowns,  CT  ,  ,  and  C,^. 

Solving  these  equations  for  the  structure  parameters,  and  applying  relations  (2) 
and  (3),  would  lead  to  the  fluxes  of  sensible  and  latent  heat. 

Of  course,  in  practice  it  is  not  as  simple  as  that.  First,  there  are  in¬ 
strumental  limitations.  These  are  discussed  in  Section  2.  A  second  limitation 
is  set  by  the  ratio  C^/C^.  If  this  ratio  is  several  times  larger  than  the  value 
assumed  in  the  above  example,  the  temperature  effect  will  dominate  at  both  0.94 
and  10.7  pm,  thus  making  the  set  of  equations  non-unique.  This  is  likely  to  occur 


Table  1.- 


-2  -3  2 

-Scintillation  parameters  at  <T>  =  298  K,  <Q>  =10  kg  m  ,  C„ 

_  ,,,  v2  -2/3  _  2  _  ,,7  . _-6  ,  2  -6  -2/3 

0.367  Km  ,  Cn  =  0.367x10  kg  m  m 


A  (pm) 

Effect  of 

temperature 

Effect  of 

humidity 

Effect  of 

temperature- 

humidity 

Transmission 

200  m 

aq 

0.94 

286 

1 

34 

0.90* 

-2.62xl0_4 

-5.79x10' 

10.7 

22 

1 

9.4 

0.96 

-2 . 16x1 0~4 

-1.55x10' 

337 

0.016 

1 

-0.26 

0.05 

-2 . 84xl0~4 

7.43x10' 

*  From  an  unpublished  calculation  of  Lawrence  and  Ochs. 


if  the  evaporation  rate  is  severely  limited  by  the  lack  of  soil  moisture,  a  situa¬ 
tion  often  encountered  in  the  non- irrigated  parts  of  eastern  Colorado. 

3 

From  (2)  and  (3),  one  immediately  infers  that  C,.,/C.  =  H/LE  x  L/p  c  ;  2.5  x  10  B, 

^  2  °2  ™ 

where  B  =  H/LE  is  the  Bowen  ratio.  The  relation  between  /C,,  and  B  has  been 
o  l  0 

pointed  out  earlier  by  Wesely  (1976).  A  Bowen  ratio  of  0.5,  like  the  one  assumed 
in  the  example,  is  characteristic  of  a  freely  evaporating  vegetation.  If  the  water 
supply  is  restricted,  (3  >  1  is  more  likely.  If  the  ratio  C,^,/C^  is  very  small,  as  may 
occur  above  a  water  surface,  humidity  fluctuations  may  be  dominant  in  the  far 
infrared,  and  even  contribute  to  the  scintillation  at  visible  wavelengths.  Finally, 
the  correlation  between  temperature  and  humidity  fluctuations  adds  to  the  variety 
of  possible  cases.  This  has  been  discussed  in  detail  by  Wesely  (1976).  A  strong 
effect  of  humidity  fluctuations  on  the  refractive  index  was  observed  by  Friehe 
et  al.  (1975). 


2.  SCINTILLATION  MEASUREMENTS— INSTRUMENTS  AND  RESULTS 

2.1  General  Features  of  the  Three  Scintillometers 

Instrumental  parameters  of  the  three  scintillometers  we  used  are  summarized 
in  Table  2.  Of  the  available  suitable  infrared  light  sources,  we  decided  to  use 
a  C02  laser  (10.7  pm)  and  an  HCN  laser  (337  pm).  The  0.94  pm  scintillometer  has 
a  LED  light  source.  As  discussed  in  the  foregoing,  we  expected  that  each  of  the 
three  contributions  to  the  refractive-index  fluctuations  (temperature  fluctuations, 

Table  2.  —  Instrument  parameters  of  the  three  scintillometers. 


A  (pm) 

Diameter 

transmitter 

(m) 

Diameter 

receiver 

(m) 

Output 

power 

(mW) 

Modulation 

frequency 

(Hz) 

Divergence 

transmitter 

(mrad) 

Divergence 
receiver 
(mrad ) 

Fresnel 

zone 

(m) 

0.94 

0.05 

0.05 

20 

25000 

6 

* 

0.0137 

10.  7 

0.002 

0.05 

3000 

260 

10 

8 

0.046 

337 

0.03 

0.25 

30 

150 

JO 

12 

0.26 

*  The  detector  is  positioned  in  the  exit  pupil  of  a  binocular  optics. 
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humidity  fluctuations,  and  correlated  temperature-humidity  fluctuations)  should  be 
prominent  in  one  of  the  scintillometers. 

The  scintillometers  were  operated  on  Table  Mountain  over  a  path  of  200  m  at  a 
height  of  4  m.  The  path  length  was  the  result  of  a  trade-off  between  the  absorp¬ 
tion  of  the  337  pm  radiation  (see  Table  1)  and  the  strength  of  the  scintillation 
(increasing  with  path  length,  see  Eq .  (5)).  The  height  was  such  that  no  outer- 

scale  effects  were  anticipated  (e.g..  Hill  and  Ochs,  1978).  Since  the  scintillation 
2 

a  decreases  with  height  (Wyngaard  and  Clifford,  1978),  it  is  wise  to  choose  the 

X 

path  height  as  low  as  possible,  giving  consideration  to  outer-scale  effects.  Table 
Mountain  is  a  flat  and  uniform  terrain  about  10  miles  north  of  Boulder.  The 
"mountain"  is  about  50  m  higher  than  its  surroundings  and  is  covered  with  a  sparse 
vegetation  of  grass.  The  scintillometers  were  set  up  near  the  optical  building, 
in  an  approximately  N-S  orientation  with  the  detectors  south,  about  10  m  from 
the  building,  on  top  of  a  tripod  (Fig.  1).  The  transmitters  were  housed  inside 
and  on  top  of  a  trailer.  In  Fig.  2  the  position  of  the  scintillometers  and  that 
of  o“her  instruments  is  indicated. 

The  output  of  the  transmitters  was  modulated  for  better  discrimination  against 
low-frequency  background  fluctuations;  besides,  the  pyroelectric  detectors  used 
for  the  C02  and  HCN  scintillometer  are  only  sensitive  to  irradiation  changes. 

The  modulation  frequencies  of  260  Hz  (10.7  pm)  and  150  Hz  (337  pm)  were  chosen 
on  considerations  of  the  high-frequency  roll-off  of  the  detectors  signal  bandwidth, 
mutual  interference,  and  60-Hz  interference.  The  0.94  pm  LED  scintillometer  can  be 
operated  on  a  much  higher  modulation  frequency  (25  kHz)  because  of  different  de¬ 
tector  characteristics.  Each  detector  unit  has  two  circular  apertures  and  two 
detector  cells.  The  receiver  apertures  of  the  C02  and  HCN  scintillometers  were 
tangent;  those  of  the  LED  scintillometer  were  not.  Besides  the  advantage  of  a 
better  rejection  of  large  eddy  sizes  (Wang  et  al.,  1978),  taking  the  difference 
of  signals  from  two  receiving  apertures  provides  compensation  for  60-Hz  pickup,  and 
for  light  source  and  background  fluctuations.  For  the  C09  and  the  HCN  scintillom¬ 
eters  pyroelectric  detectors  were  used  because  they  are  easy  to  handle  (no  cooling 
required)  and  provide  sufficient  sensitivity  for  our  purpose  (at  least,  in  theory). 

As  the  arrangement  of  the  scintillometers  differs  from  a  point  source/point 
detector  setup,  the  expression  relating  the  scintillation  to  the  structure  param¬ 
eters  differs  from  (5)  by  a  factor  depending  on  the  geometry  of  the  apertures. 
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Figure  1. — Scintillation  detectors.  The  lower  box  is  the  HCN  receiver;  it  con¬ 
tains  two  25-cm  f/1  parabolic  mirrors.  The  smaller  box,  mounted  to  the  left 
of  the  lower  box,  contains  the  electronics.  On  top  of  the  HCN  receiver  are 
the  CO^  receiver  (left)  and  the  LTD  receiver  (right). 


Figure  2. — Section  of  Table  Mountain,  indicating  the  positions  of  the  transmitters 
(T),  receivers  (R) ,  the  trusts  with  the  in  situ  instruments  (M)  and  the  optical 
building  (OB)  where  the  signals  were  processed. 
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Furthermore,  it  is  necessary  to  add  a  correction  factor  that  accounts  for  the  contri¬ 
bution  of  eudies  near  the  inner  scale  of  the  turbulence,  because  the  refractive- 
index  spectrum  exhibits  a  bump  at  the  corresponding  spatial  wave  numbers  (Hill, 
1978a).  Taking  these  factors  into  account,  the  expressions  for  the  scintillation 
are  these: 


2  i  a  o  9  T 

0.94  pm:  a  =  1.54  x  10  C  (0.809  +  0.920  £  -  0.161  £  -  0.177  £  ) 

n  o  o  o 

(0.4  cm  <  £  <  1.1  cm) 

~  0  - 

2  1  n  9  9 

10.7  pm:  a  =  2.830  x  10  C  (0.824  +  1.279  £  -  0.736  £  ) 

n  o  o 

(0. 3  cm  <  £  <1.1  cm) 

~  o  ~ 

337  pm:  a2  =  5.380  x  108  C  2  (1.011  -  0.054  £  +  0.088  £  2) 

n  o  o 

(£q  <1.1  cm) 


(9) 


(10) 


(ID 


where 


2  <(IA  -  V  > 
°  ^ - 


(12) 


2  _  -2/3 

Here  is  in  units  m  ;  £q  is  the  inner  scale  of  turbulence  in  centimeters; 

£f  =  7.4  n.  where  n  is  the  Kolmogorov  microscale.  (See,  e.g.,  Hill  and  Clifford, 
1978.)  1^  and  lfi  are  the  light  intensities  observed  in  circular  apertures  labeled 

A  and  B.  It  is  assumed  that  <I^>  =  <-*-g>  =  <T>.  The  apertures  of  the  CO^  and  HCN 
transmitters  are  assumed  to  be  infinitesimally  small  (point  sources);  the  LED 
transmitter  aperture-averaging  effect  must  be  taken  into  account.  Transmitter 
and  receiver  aperture-averaging  calculations  give  Eq.  (9)  for  the  LED  scintillom¬ 
eter.  Receiver  aperture  averaging  calculations  give  Eqs.  (10)  and  (11)  for  the 
CO2  and  HCN  scintillometers.  The  calculations  were  done  by  R.J.  Hill  (private 
communication).  The  expression  between  parentheses  is  the  inner-scale  correction. 
It  is  small  for  the  HCN  scintillometer  (2%  at  £q  =  0.75  cm),  but  quite  considerable 
for  the  other  two  instruments  (33%  for  the  LED,  and  37%  for  the  CO^  scintillometer 
at  ^0  =  cm).  The  inner-scale  correction  was  calculated  using  the  spectrum- 

with-bump  for  temperature  only.  This  might  not  be  accurate  enough  for  the  CO^ 
scintillometer,  where  the  correction  is  large  and  the  differences  between  the 
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o 

Figure  3. — Cn^  path-weighting  functions.  The  0.94-ym  curve  is  for  a  5-cm  tangent 
aperture  receiver,  5-cm  transmitter  scintillometer.  The  inner  scale  is 
0.75  cm  for  all  cases.  From  R.J.  Hill  (unpublished). 

temperature,  humidity,  and  temperature-humidity  spectra  may  be  of  importance. 

Caution  has  to  be  used  if  the  third  term  in  Eq.  (6)  is  negative  and  tends  to  cancel 
the  first  two  terras.  The  path-weighting  functions  for  =  0.75  cm  are  shown  in 
Fig.  3  (also  calculated  by  R.J.  Hill).  Their  shapes  are  not  very  sensitive  to 

2 

To  illustrate  the  kind  of  signal  we  expect,  0  has  been  calculated  for  the 

f  n  2  „  v2  -2/3  „  2  .  ,,,  -6  .  2  -6  -2/3  _  r~  2r  2,1/2 

case  of  CT  =  0.367  Km  ,  =  0.367  x  10  kg  m  m  ,  CTq  =  [CT  ]  , 

<T>  =  298  K  and  <Q>  =  10  ^  kg  m  The  results  are  shown  in  Table  3,  where  the 

2  1/2 

scintillation  signals  < (1^  -  Ig)  >  are  also  given.  These  signals  should  be 
compared  with  the  noise  power.  Noise  power  was  derived  from  the  manufacturer's 
data  on  the  detectors,  assuming  a  bandwidth  of  100  Hz  for  the  CO^  and  HCN  detectors. 
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Table  3. — Expected  performance  of  the  scintillometers,  based  on  Tables  1  and  2, 


and  on  formulas  (9)- (11)  with  l  =0.8  cm. 

o 


A  (pm) 

Power  on 
detector 
(W) 

<<ia"V2> 

=  s 

(W) 

Noise 

detector=N 

(W) 

S/N 

C  2  x  innei 
n 

scale  corr, 
<»-2/3) 

<I2> 

0.94 

3. lxlO-5 

74xl0-3 

8. 5xl0-6 

_ 

* 

-13 

3.95x10 

10.7 

1.8x10" 3 

ll.lxlO"3 

2. OxlO-5 

3xl0~8 

667 

-13 

3.93x10 

-5 

-3 

-7 

-8 

-11 

337 

2.4x10 

8.4x10 

2.2x10 

4x10 

5.5 

1.57x10 

*  The 

-  2 

noise  equivalent  C  of  the 

0.94  pm  LED  scintillometer  is 

lower 

in-16  -2/3 
than  10  m 

11 

It  is  seen  that  the  HCN  scintillometer  has  the  lowest  signal-to-noise  ratio. 

2 

The  state  of  affairs  represented  in  Table  3  is  a  rather  optimistic  one:  is 

very  high  in  view  of  the  often  desiccated  terrain  of  Table  Mountain,  and  no  account 
has  been  taken  of  other  sources  of  noise,  such  as  electronic  and  acoustic  noise. 


2.2  The  0.94-pm  LED  Scintillometer 

For  optical  and  mechanical  details  of  the  instrument  the  reader  is  referred 

to  a  report  by  Ochs  et  al.  (1980).  The  transmitter  was  mounted  on  top  of  the 

trailer  in  which  the  CO^  and  HCN  laser  were  located.  The  receiver  was  put  on 

top  of  the  HCN  receiver  box.  Although  initially  a  scintillometer  with  15-cm 

apertures  was  used,  it  was  replaced  by  one  having  5-cm  apertures,  after  some 

2 

preliminary  measurements.  The  reason  was  that  the  observed  C^  was  several 
times  larger  than  calculated  from  C ^  measurements  with  fast  temperature  sensors, 
and  there  was  little  experience  with  that  scintillometer.  But  also  the  5-cm 
aperture  scintillometer  gave  too  high  a  reading.  It  was  next  found  that  the 
effective  aperture  of  the  receivers  was  4  cm  rather  than  5  cm,  and  that  the  trans¬ 
mitter  aperture  was  partially  coherent.  The  first  effect  accounts  for  an  over- 
2  7 /6 

estimate  of  C  by  a  factor  of  (5/4)  '  =  1.24;  the  second  effect  made  the  reading 

n 

1.62  times  too  high.  (This  was  found  by  inserting  a  ground  glass  in  front  of  the 
LED.)  A  third  correction  because  of  the  receiver  aperture  separation  of  6.7  cm 
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r 


(center  to  center)  is  included  in  Eq.  (9).  The  reading  of  the  scintillometer  was 

2 

corrected  for  these  effects  and  the  comparison  with  the  in  situ  observed  is 

shown  in  Fig.  4.  The  correspondence  at  low  wind  speed  is  quite  satisfactory,  but 

at  higher  wind  speeds  (4  ra  s  ^  has  been  used  as  the  criterion)  the  optically  ob- 
2 

served  is  too  high.  Most  likely,  this  is  caused  by  vibrations  of  the  detector 

box  at  higher  wind  speeds. 


2.3  The  10.7-ym  CO^  Scintillometer 

A  Sylvania  model  950  CO^  laser  served  as  radiation  source.  The  laser  was 
tuned  to  the  P-30  line  (10.697  ym)  and  had  a  continuous-wave  output  power  of  5  W. 
Part  of  the  light  beam  was  reflected  by  a  ZnSe  plate  to  a  thermopile  radiation 
detector  for  monitoring  the  output  power  of  the  laser.  Then  the  beam  was  com¬ 
pressed  by  means  of  two  germanium  lenses  to  one  third  of  its  original  diameter, 
thus  increasing  the  beam  divergence  by  a  factor  of  3  to  about  10  mrad  and  making 
the  aiming  of  the  beam  at  the  receiver  less  critical.  After  passing  an  attenuator 
and  a  260-Hz  mechanical  chopper,  the  beam  was  reflected  upwards,  led  through  an 
aluminum  tube  through  the  roof  of  the  trailer,  and  directed  to  the  receiver  by 
means  of  a  second  mirror  mounted  on  an  adjustable  support.  Both  mirrors  were 
aluminum,  front-surface  plane  mirrors.  The  distance  from  the  beam  compressor  to 
the  second  mirror,  as  measured  along  the  light  path,  was  2  m. 

The  receiver  for  the  10. 7-ym  scintillometer  was  mounted  on  top  of  the  receiver 
for  the  337-ym  scintillometer,  next  to  the  0.94-ym  receiver.  The  two  tangent 
circular  apertures  each  have  a  diameter  of  5  cm  (which  is  nearly  equal  to  the 
Fresnel  zone  size),  and  were  covered  with  10-ym-thick  plastic  foil  to  prevent  air 
currents  from  entering  the  box.  It  was  found  that  air  currents  seriously  degrade 
the  performance  of  the  detector  cells.  The  radiation  was  focused  by  two  spherical 
mirrors  (f  =  28  cm)  onto  two  pyroelectric  detectors  (Molectron  Pl-72)  having  a 
circular  sensitive  area  of  2  mm  diameter.  The  detectors  were  protected  by  germanium 
windows  with  an  anti-reflection  coating.  The  output  signal  of  each  detector  was 
converted  to  a  low  impedance  by  a  unity-gain  amplifier  mounted  near  the  detector, 
and  was  then  led  through  a  two-stage  amplifier  (maximum  gain  10,000  x)  and  a  band¬ 
pass  filter  (gain  1.8  x,  3-dB  points  at  150  and  400  Hz),  which  were  situated  in  a 
box  mounted  directly  under  the  box  containing  the  mirrors  and  detectors.  The  final 
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Figure  4. — Cnz  from  the  0.94-ym  LED  scintillometer,  compared  with  Cn2  calculated  from  CT2  measurements. 
The  dots  are  for  windspeeds  lower  than  4  m  s”l,  the  crosses  for  higher  windspeeds. 


output  signal  was  maintained  at  a  constant  rms  value  by  means  of  an  automatic  gain 
control  feedback  loop  having  a  time  constant  of  about  Is.  A  schematic  of  the 
electronics  is  given  in  Appendix  1.  The  electronic  filters  were  chosen  so  that  the 
amplitude-modulated  260-Hz  carrier  wave  would  pass  freely,  while  other  signals 
contributing  only  to  the  noise  of  the  system  were  suppressed.  The  scintillation 
signal  is  expected  to  have  dominant  contributions  from  aperture-sized  eddies,  which 
are  convected  across  the  scintillation  path  by  the  wind.  Thus,  at  a  wind  speed  of 
5  m  s  \  the  frequencies  of  the  scintillation  are  centered  around  5ms  ^/5  cm  = 

100  Hz. 


The  two  output  signals  were  processed  by  a  Nova  computer  inside  the  building. 
During  16  minutes,  600  samples  of  each  signal  were  taken,  say  A  and  B.  Then  the 
following  quantity  was  calculated: 


S 


A 

<(|A|  -  ] B  j  < 

B 

>}  > 

<A2> 


(13) 


By  this  procedure,  the  signals  were  demodulated  and  corrected  for  any  difference 

in  mean  value  (by  the  factor  < | A | >/< | B | >) .  The  angle  brackets  indicate  the 

2 

averages  of  the  600  samples.  S  is  the  same  quantity  as  a  in  Eq.  (12).  This 
procedure  was  followed  instead  of  an  analog  computation  (as  was  done  with  the 
0.95-pm  scintillometer),  because  the  analog  demodulation  of  the  signal  offered 
problems  caused  by  the  large  ratio  of  bandwidth  to  carrier  frequency. 


The  computer  has  a  time  lag  of  40  us  between  the  sampling  of  two  neighboring 

channels.  This  was  compensated  by  an  approximately  equivalent  phase  shift  applied 

to  one  of  the  input  signals,  by  means  of  an  RC  filter.  The  desired  phase  shift 

was  0  =  -  40  x  10  ^  u)  while  that  provided  by  the  RC  filter  (RC  =  40  x  10  ^  sec) 

was  ip  =  Arctan  (-  40  x  10  ^  ..j)  .  The  largest  u.  desired  is  very  small  compared  to 
— 6  __ 

(40  x  10  )  rad/sec;  thus  the  two  phase  shifts  are  nearly  equal  because  the 

arctangent  is  approximately  linear  in  its  argument  provided  the  argument  is  much 
smaller  than  unity.  The  proper  functioning  of  the  procedure  was  checked  with 
modulated  sine-wave  signals  from  sine-wave  generators.  With  no  modulation  at  all, 
the  value  of  S  was  =  2  x  10^,  a  value  thought  to  be  caused  by  the  10-raV  noise  of 
the  a-d  converter.  The  noise  level  of  S  was  sufficiently  low  compared  with  the 

_3 

expected  signal  level  (10  ). 
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Measurements  of  a  (=  S)  were  compared  with  values  calculated  from  in  situ 
2  2 

measurements  of  ,  and  C^q,  using  fast  temperature  and  humidity  sensors. 

The  contribution  of  C  2  to  C  ^  was  negligible,  and  that  of  C_,,  at  most  20%.  The 

Q  n  TQ  2 

results  are  shown  in  Fig.  5.  It  is  found  that  the  observed  values  of  o  are 
larger  than  the  calculated  values  by  a  factor  1.5-2.  Part  of  this  discrepancy 
may  be  ascribed  to  instability  of  the  pyroelectric  detectors.  It  was  noticed 
that  the  amplitude  of  the  received  signal  showed,  from  time  to  time,  irregular 
changes  that  were  occasioned  by  neither  the  transmitter  nor  the  atmosphere.  The 


TQ' 
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changes  were  much  faster  than  the  response  time  of  the  automatic  gain  control 
loop.  It  is  thought  that  a  thermal  imbalance  between  the  detector  cell  and  its 
environment  might  be  responsible  for  this  behavior,  but  time  was  lacking  to  in¬ 
vestigate  this  problem.  As  a  result,  few  conclusions  can  be  drawn  from  the  mea¬ 
surements  with  the  CO ^  scintillometer. 


2.4  The  337-ym  Far-Infrared  Scintillometer 

The  HCN  laser  that  was  used  as  a  radiation  source  has  a  4-m-long  cavity  with 
a  folded  confocal  geometry  (one  flat  mirror,  and  one  concave  mirror  with  a  focal 
length  of  8m).  As  an  output  coupler,  a  polyethylene  beam  splitter  arranged  in 
a  Michelson  interferometer  configuration  was  used.  Output  coupling  could  thereby 
be  varied  between  0%  and  16%.  Details  of  the  laser  can  be  found  in  a  report  by 
Wells  et  al.  (1971).  Because  of  the  cavity  geometry,  the  output  beam  should  have 
a  diameter  of  4.1  cm  close  to  the  laser  and  finally  reach  a  full-angle  divergence 
of  10  mrad. 

The  laser  was  powered  by  a  regulated  3Q00-V  d.c.  power  supply  with  a  ballast 

resistor  of  1  kft  in  series.  At  a  discharge  current  between  0.7  and  1  A,  the 

3 

anode-cathode  voltage  is  ^2000  V.  A  constant  flow  of  20  cm  /min  (STP)  of  a  mixture 
of  He,  N^,  and  CH^  (in  the  ratio  3:36:8)  is  maintained  through  the  Pyrex  discharge 
tube.  The  lasing  HCN  gas  is  produced  by  chemical  reactions  in  the  discharge.  A 
constant  gas  pressure  in  the  tube  is  rather  easily  produced  by  equilibrium  between 
the  gas  supply  rate  (provided  by  a  gas  cylinder  filled  with  the  proper  mixture,  and 
a  fine  needle  valve)  and  the  pumping  rate  (provided  by  a  mechanical  vacuum  pump). 

The  composition  of  the  gas  fed  to  the  laser  was  chosen  on  the  basis  of  a 
publication  of  Belland  and  Veron  (1973),  who  claimed  that  when  this  mixture  was 
used  the  discharge  tube  remained  clear,  in  contrast  to  the  brown  deposit  that 
usually  forms  on  the  wall,  and  the  output  power  was  high.  While  our  discharge 
tube  stayed  clean,  the  output  power  was  less  than  expected  on  the  basis  of  calcu¬ 
lations  and  experiments  of  Belland  et  al.  (1976).  This  may  be  attributed  to  the 
non-optimum  tube  diameter  (8  cm),  and  non-optimum  electrode  geometries.  The 
highest  output  power,  obtained  at  a  gas  pressure  of  1.2  torr,  was  ^50  mW.  The 
discharge  was  then  rather  turbulent,  however.  At  even  higher  gas  pressures  the 
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power  supply  could  not  meet  the  higher  voltage  demand  of  the  discharge;  consequently, 
the  current  dropped,  and  so  the  lasing  power.  At  lower  gas  pressures  the  discharge 
became  less  turbulent,  until  at  “M).9  torr  a  stable,  striated  discharge  was  obtained 
with  ^30  mW  in  the  337  pm  line.  (Line  identification  was  simple  with  the  Michelson 
interferometer  output  coupler.)  After  a  considerable  warming-up  period  (^1  hour) 
the  laser  usually  ran  very  stably. 

An  effort  was  made  to  get  a  stable  discharge  with  a  higher  laser  power.  Other 
gas  mixtures  were  tried  (our  mix  with  CH^  or  NH^  added,  and  a  mix  of  CH^  and  NH^ 
only),  the  gasflow,  discharge  current,  and  gas  pressure  were  varied,  and  other 
electrode  geometries  were  used,  but  all  with  little  or  no  improvement. 

An  effort  was  made  to  measure  the  beam  diameter  and  divergence.  To  do  this, 
apertures  were  put  in  the  beam  at  positions  between  0.13  and  10.1  m  from  the  exit 
window.  Directly  after  the  aperture,  the  radiation  was  focused  by  a  polyethylene 
lens  on  a  thermopile  detector.  In  this  way  the  beam  diameter  was  sensed.  It  was 
observed  that  the  beam  diameter  at  the  exit  window  was  ^2  cm,  instead  of  the  4  cm 
calculated,  and  that  the  beam  divergence  was  about  the  expected  value,  10  mrad. 
However,  about  a  50%  radiation  loss  was  observed  at  10.1  m,  which  cannot  be  ex¬ 
plained  by  divergence  nor  by  absorption. 

The  laser  beam  was  chopped  at  150  Hz,  led  through  an  aluminum  tube  through 
the  roof  of  the  trailer  and  aimed  at  the  receiver  in  much  the  same  way  as  for  the 
CO 2  laser  beam.  A  plastic  beam  splitter  reflected  part  of  the  beam  onto  a  Teflon 
lens,  which  focused  the  radiation  on  a  pyroelectric  detector  that  served  as  a 
monitor.  The  distance  from  the  laser  exit  window  to  the  mirror  on  top  of  the 
trailer  was  3. 35  m. 

The  receiver  unit  was  a  box  containing  two  25-cm-diameter  tangent  mirrors  and 
two  pyroelectric  detectors.  The  mirrors  had  a  focal  length  of  25  cm,  and  were 
machined  from  cast  aluminum  to  a  parabolic  shape  by  the  NBS  mechanical  workshop. 

The  pyroelectric  detectors  (Molectron  Pl-73)  had  a  sensitive  circular  area  of  3  mm 
diameter  and  Teflon  windows  (0.8  mm  in  thickness).  The  radiation  entered  the 
detector  box  through  another  Teflon  window,  1.5  mm  thick.  By  these  windows,  an 
effective  blocking  of  the  CO2  radiation  was  achieved:  the  transmission  at  10.7  pm 
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was  estimated  at  10  at  337  pm  it  was  0.7.  The  area  of  the  pyroelectric  de¬ 
tectors  was  large  enough  to  accommodate  for  mirror  imperfections  (the  image  of  a 
point  source  is  blurred,  and  has  a  diameter  of  M  mm),  and  for  a  possible  tilt  of 
the  wavefront  at  the  receiver  (which  might  let  the  focal  spot  wander  over  at  most  1 
mm)  . 

Signal  processing  was  along  the  same  lines  as  for  the  CO^  scintillometer. 

The  signals  were  first  amplified  (gain  2x)  near  the  detectors,  then  led  to  a  small 
box  mounted  to  the  box  containing  the  mirrors  and  detectors,  where  a  maximum  gain 
of  2.5  x  10^  was  applied;  finally  the  signals  passed  a  bandpass  filter  (3-dB  points 
at  80  and  240  Hz).  The  total  maximum  gain  was  1.7  x  10"*,  and  the  final  output 
signals  were  maintained  to  a  constant  root-mean-square  value  by  means  of  an  auto¬ 
matic  gain  control.  The  expected  maximum  bandwidth  of  the  signal  was  40  Hz.  In 
order  to  suppress  60-Hz  pickup,  the  electronics  were  powered  by  batteries  placed 
inside  the  box.  The  electrical  schematic  is  given  in  Appendix  I . 

In  operating  the  HCN  scintillometer,  it  was  found  that  the  system  was  severely 

hampered  by  instability  of  one  of  the  pyroelectric  detectors.  (The  same  problem 

was  found  with  the  CO^  scintillometer,  but  to  a  lesser  extent.)  Most  of  the 

time,  no  measurement  could  be  made  at  all.  On  only  one  day  did  the  system  operate 

-4  2 

reasonably  well  for  a  couple  of  hours.  A  value  of  7  x  10  was  observed  for  0  ; 

the  calculated  value  was  8  x  10  b.  Obviously,  only  the  noire  of  the  system  at 

2 

that  particular  moment  of  time  was  observed.  The  low  calculated  value  of  o 

indicates  the  kind  of  problem  one  faces  in  trying  to  observe  scintillation  with 

submillimeter  waves  in  a  very  dry  atmosphere:  it  is  often  too  low  to  be  detectable. 

During  a  period  of  three  months  (January  through  March  1981)  on  only  one  day 

(2  February)  were  the  humidity  fluctuations  strong  enough  to  give  an  observable 
2  -3 

j  =  1.2  x  10  .  Unfortunately,  on  that  day  the  scintillometer  did  not  work. 

It  is  suggested  that  before  making  new  at t^^pto  to  observe  scintillation 
at  337  pm,  one  should  first  have  a  thorough  look  at  the  functioning  of  pyroelectric 
detectors  in  a  poorly  temperature-controlled  environment.  If  the  use  of  such  a 
detector  is  too  marginal,  other  devices  should  be  considered,  like  a  germanium 
bolometer,  though  these  detectors  have  to  be  cooled  to  cryogenic  temperatures, 
which  offers  some  technical  complications. 


3.  IN  SITU  MEASUREMENTS — THE  INSTRUMENTS 


These  micrometeorological  in  situ  measurements  were  acquired  to  serve  as  a 
reference  for  the  optical  measurements,  and  to  give  additional  meteorological 
information: 


(a)  Wind  speed  and  wind  direction  with  a  Gill  propeller  vane. 

(b)  Air  temperature  and  dewpoint  temperature  with  an  EG&G  model  220  system. 

(c)  Net  radiation  with  a  Swissteco  radiometer. 

(d)  Temperature  fluctuation,  and  temperature  difference,  using  three 
2-um-diameter  Pt  wires. 

(e)  Humidity  fluctuations,  using  two  Ly-a  hygrometers. 

(f)  Vertical  wind  speed  fluctuations,  with  a  20-cm-path  sonic  anemometer. 

An  example  of  the  arrangement  of  the  sensors  (a),  (d),  (e),  and  (f)  is  given 
in  Fig.  6.  These  sensors  were  mounted  at  a  height  of  4  m  on  top  of  two  well-guyed 
poles  situated  at  the  midpoint  of  the  optical  path.  The  sensors  were  located  so 
that  the  following  quantities  could  be  measured: 

2 

(a)  CT  ,  from  the  temperature  difference  of  two  Pt  sensors. 

2 

(b)  Cq  ,  from  the  humidity  difference  of  two  Ly-a  hygrometers. 

(c)  from  one  pair  of  temperature  sensors  and  one  pair  of  humidity  sensors. 

(d)  Sensible  heat  flux,  by  correlating  the  vertical  wind  speed  with  the 
temperature  fluctuation  measurement. 

(e)  Latent  heat  flux  (or  moisture  flux),  by  correlating  the  vertical  wind 
speed  with  humidity  fluctuations. 


The  orientation  of  the  separation  between  the  two  Ly-a  hygrometers  was  manually 
set  to  near  perpendicularity  to  the  wind  direction.  That  means  that  cross-wind 
structure  parameters  were  observed.  There  is  evidence  that  there  is  no  difference 
between  crosswind-  and  streamwise-measured  structure  parameters  (Wyngaard  et  al., 
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1978).  No  experimental  evidence  is  available  regarding  vertically  measured  struc¬ 
ture  parameters. 


The  distance  between  the  Pt  wires  and  the  center  oi  the  Ly-  i  humidiometers 
was  set  at  3.7  cm.  The  wires  were  positioned  upwind  from  the  Ly-i's.  It  is 
assumed  that  the  presence  of  the  l.y-u's  does  not  significantly  disturb  the  adverted 
temperature  field.  Some  support  of  this  assumption  was  found  by  observing  the 
temperature  variance  at  various  distances  between  a  temperature  sensor  and  a  i.y-  i . 
Another  Pt  wire  was  at  a  fixed  position  and  served  as  a  reference.  Over  the 
whole  range  of  available  distances,  10.5  to  2.5  cm,  no  difference  was  noted.  in 
another  experiment,  two  temperature  spectra  were  measured  sequentially,  one  with 
the  sensor  at  4  cm  from  the  Ly-  .  axis,  the  other  with  the  sersor  on  the  axis,  and 
thus  in  the  gap  between  the  source  and  detector.  The  spectra  were  similar. 

The  distance  between  the  axis  of  the  sonic  anemometer  and  the  temperature 
and  humidity  sensor  to  be  correlated  with  vertical  wind  was  10  cm  at  most. 

The  propeller  vane  was  situated  close  to  the  turbulence  instruments  because 
the  windspeed  is  needed  for  subsequent  data  analysis,  where  Taylor’s  hypothesis 
is  used  for  the  calculation  of  structure  parameters  from  time-lagged  measurements 
(see  Section  4).  The  other  instruments — the  net  radiometer  and  the  temperature/ 
dewpoint  instrument — were  situated  at  some  distance  from  the  poles.  The  tempera¬ 
ture,  dewpoint  temperature,  and  net  radiation  measurements  served  for  obtaining 
general  meteorological  information.  The  dewpoint  temperature  measurements  were 
further  used  for  calibrating  the  Ly-o  humidiometers,  as  is  discussed  below.  All 
instruments  were  set  up  permanently  in  the  field,  with  the  exception  of  the  Ly- . 
hygrometers  and  the  Pt  temperature  sensors.  The  signals  were  amplified  close  to 
the  instruments  before  being  transported  to  the  Nova  computer  inside  the  building 
(the  only  exception  being  the  signal  of  the  net  radiometer,  which  was  directly 
connected  to  the  building).  There,  a  final  amplification  was  performed  to  ac¬ 
commodate  the  signal  level  to  the  +  5-V  input  range  of  the  computer.  Special  care 
had  to  be  taken  to  avoid  ground  loops  between  the  computer  ground  and  the  instru¬ 
ments. 


Next,  we 
for  using  the 


shall  pay 
Ly-  -  line 


some  special  attention  to  the  i.y- 
(121.6  nni)  for  measuring  humidity 


hygrometer.  The  reason 
is  that  the  ratio  of 


water  vapor  and  oxygen  absorption  shows  a  remarkable  peak  very  close  to  the  Ly-u 
line  (Tillman,  1965).  Furthermore,  the  absorption  is  fairly  strong,  thus  making 
measuring  over  a  short  path  (M  cm)  possible.  Another  important  consideration  is 
that  the  construction  of  a  Ly-u  source,  a  hydrogen  discharge  tube,  is  relatively 
easy.  A  nitric  oxide  ionization  chamber  is  usually  used  as  a  detector.  Both 
source  and  detector  should  have  special  windows  for  transmitting  the  Ly-a  radia¬ 
tion.  The  best  choice  seems  to  be  magnesium  fluoride,  which  resists  the  weather 
(especially  moisture-resistant)  fairly  well.  Magnesium  fluoride  is  opaque  below 
115  nm,  while  the  ionization  potential  of  nitric  oxide  is  at  135  nm.  The  de¬ 
tector  is  thereby  sensitive  to  the  wavelength  interval  115-135  nm  only,  which 
nicely  brackets  the  Ly-u  line  (121.6  nm).  Since  no  other  means  of  spectral  filter¬ 
ing  can  be  used  in  the  field,  it  is  important  to  have  the  emission  of  the  Ly-u 
line  as  pure  as  possible  in  order  to  minimize  the  absorption  by  spurious  oxygen 
lines.  Spectral  purity  can  be  controlled  by  the  choice  of  gas  mixture  and  pressure 
in  the  discharge  tube.  An  important  refinement  is  the  addition  of  a  UH^  re¬ 
servoir  to  the  discharge  tube.  Uranium  hydride  provides  a  constant,  low,  hydrogen 
pressure  if  heated.  In  this  way  a  very  pure  Ly-a  source  with  a  lifetime  of  more 
than  1000  hours  is  obtained  (Buck,  1977).  A  disadvantage  of  the  UH^  source  might 
be  its  larger  size  in  comparison  with  the  "conventional",  unheated  source,  which 
makes  it  more  likely  to  obstruct  the  airflow  through  the  gap  between  source  and 
detector.  However,  little  is  known  about  this  effect.  We  preferred  to  use  the 
conventional  tube. 

The  manner  in  which  the  sources  and  detectors  were  mounted  is  shown  in  Fig.  7. 

A  small  box,  containing  two  amplifiers,  was  situated  about  0.5  m  from  the  detectors. 
(For  the  electronic  schematic,  see  Appendix  II.)  The  output  voltage  was  about  1  V. 
The  frequency  response  of  the  amplifiers  was  d.c,  to  500  Hz  (3-dB  point).  A 
laboratory  test,  where  the  frequency  response  of  the  Ly-a  hygrometer  was  measured 
by  electronic  modulation  of  the  source  current,  showed  that  these  amplifiers 
essentially  limited  the  over-all  response  of  the  hygrometer.  With  faster  ampli¬ 
fiers,  an  almost  flat  frequency  response  from  d.c.  to  10  kHz  was  found  (J.T. 
Priestley  and  W.D.  Cartwright).  The  output  voltage  fluctuations  (reflecting  hu¬ 
midity  fluctuations)  we  observed  were  often  so  small  that  it  was  necessary  to 
offset  the  output  voltage  and  amplify  it  again  in  order  to  boost  the  fluctuating 
signal  sufficiently  above  the  computer  noise  (M0  mV).  The  noise  of  the  hygrom¬ 
eter,  as  measured  at  the  output  of  the  first  amplifier,  was  less  than  1  mV  peak-to- 
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Figure  7. — Construction  of  a  Ly-a 
hygrometer. 


_3 

peak,  which  corresponds  to  a  humidity  noise  of  0.01  g  m  or  less.  The  offset  of 
the  amplifiers  can  easily  be  checked  by  blocking  the  light  path.  It  should  be  less 
than  1  mV.  We  found  that  it  had  to  be  adjusted  only  very  rarely. 


The  humidity  was  calculated  from  the  output  voltage  of  the  hygrometers,  using 
a  procedure  of  Buck  (1980).  The  procedure  is  based  on  laboratory  calibrations  of 
similar  sources  and  detectors,  and  accounts  for  additional  oxygen  absorption. 

Also,  a  correction  factor  for  the  changing  collimation  with  different  distances 
between  source  and  detector  is  included.  The  calculation,  which  was  performed  on 
the  Nova  computer,  also  involved  a  calibration  "constant".  This  constant  was 
obtained  by  comparison  with  an  absolute  hygrometer,  in  our  case  the  dewpoint  in¬ 
strument.  It  was  necessary  to  adjust  this  calibration  constant  frequently  because 
drift  of  the  Ly-a  hygrometers  could  be  as  high  as  10%  per  hour.  To  do  so,  each 
computer  sample  of  the  two  Ly-a  hygrometers  was  compared  with  the  reading  of  the 
dewpoint  hygrometer,  and  the  constant  was  adjusted  slightly  (fractional  change  5  x 
10  “*).  By  this  procedure  a  drift  rate  of  10%  per  hour  could  be  compensated,  with¬ 
out  affecting  the  fluctuations. 


The  limited  experience  we  have  gained  with  the  Ly-a  hygrometer  as  an  outdoor 
instrument  is  not  positive  in  all  respects.  Although  the  instrument  should  be 
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rugged  enough  to  withstand  all  kinds  of  weather,  it  was  found  that  the  lifetime 
of  the  detectors  was  less  than  the  several  hundred  hours  anticipated.  Detector 
failure  was,  in  all  cases,  occasioned  by  leaks  in  the  sealing  between  the  magnesium 
fluoride  window  and  the  Pyrex  tube.  In  one  case,  the  detectors  failed  after  an 
exposure  of  several  weeks  to  varying  weather,  including  freezing  rain  and  snow, 
with  only  20  hours  of  operation.  The  failure  was  manifested  by  an  initially  low 
detector  resistance,  followed  by  a  high  resistance  but  low  sensitivity.  After 
this  experience,  it  was  decided  to  put  the  instrument  outside  only  during  mea¬ 
surements.  Since  then,  the  hygrometers  have  performed  more  satisfactorily. 


4.  DATA  PROCESSING  AND  ON-LINE  CALCULATIONS 

Data  processing  was  done  on  a  Nova  computer.  The  computer  had  no  facility 
for  storing  data  on  disc  or  magnetic  tape.  Therefore,  the  calculation  of  mean 
values,  variances,  structure  functions,  etc. ,  had  to  be  done  on-line.  A  dis¬ 
advantage  of  this  system  is  that  it  is  not  possible  to  calculate  spectra.  For  a 
few  runs,  however,  we  recorded  data  of  the  fast  temperature  sensor  and  one  Ly-a 
hygrometer  on  an  analog  tape  recorder  and  digitized  these  data  later.  An  advantage 
of  the  on-line  calculations  is  that  there  is  an  immediate  check  on  sensor  perform¬ 
ance. 


The  computer  was  run  with  either  of  two  programs:  a  program  called  METEO, 
handling  nine  input  channels  with  a  sampling  rate  of  142  Hz,  and  a  program  called 
SCINT,  handling  fourteen  channels  at  50  Hz.  The  sampling  rates  are  for  both  cases 
the  highest  possible.  The  reason  for  having  a  high  sampling  rate  will  be  made 
clear  below.  The  input  channels  for  METEO  are  the  two  Ly-a  hygrometers,  fast 
temperature,  fast  temperature  difference,  sonic  anemometer,  net  radiation,  wind 
speed,  temperature,  and  humidity.  SCINT  handles  the  same  inputs,  plus  one  channel 
for  the  0.94  ym  scintillometer,  and  two  channels  for  each  of  the  other  two  scin¬ 
tillometers. 

The  high  sampling  rate  was  for  calculating  structure  parameters  from  time- 
lagged  observations.  For  instance,  the  observation  of  the  humidity  Q  at  time 
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points  t  and  t  +  T  enables  one  to  calculate  the  humidity  structure  parameter  by 
°  °  2  -2/3 

<[Q(tQ)  -  Q(tQ  +  t)]  >(<u>t)  ,  where  <u>  is  the  mean  windspeed  for  the  period 

of  observation  considered.  The  shortest  possible  time  delay  T  is  obviously  set 
by  the  sampling  rate  of  the  computer.  The  purpose  of  these  measurements  is  to 
investigate  the  supposed  independence  between  structure  parameters  and  distance 
r  by  simply  varying  the  time  delay.  The  alternative  method,  an  array  of  many 
sensors,  is  less  attractive. 

In  more  detail,  the  data  processing  is  as  follows.  The  computer  takes  in 
data  during  1  second.  Using  program  SCINT,  the  first  sample  of  the  142  samples  for 
each  channel  is  stored,  and  only  for  the  fast  temperature  channel  and  one  Ly-a 
channel  the  samples  3,  3+21  (i=0, 1, 2, . . . ) ,  up  to  131  and  142  are  retained.  The 
latter  samples  serve  for  calculating  time-lagged  structure  parameters.  With  program 
METEO,  the  samples  3,  4,  5,  7,  11,  19,  35,  and  50  are  stored  for  this  purpose. 

The  temperature  and  humidity  differences  are  related  to  the  third  sample  because  it 
was  found  that  at  a  sampling  rate  of  142  Hz  the  time  interval  between  the  first 
and  the  second  sample  was  longer  than  1/142  s,  and  the  interval  between  the  second 
and  the  third  sample  was  too  large  at  times.  After  all  data  are  transferred  into 
physical  quantities,  a  number  of  calculations  are  performed  (like  squaring,  cross 
products,  etc.),  and  the  results  are  stored.  This  takes  about  0.5  second.  Then, 
the  next  1-second  block  of  data  is  taken  in.  After  a  run  of  600  blocks,  which 
takes  about  16  minutes,  mean  values  are  calculated  and  printed  on  a  teletype. 

Next,  we  give  formulas  that  were  used  for  the  calculation  of  variances, 
structure  parameters,  etc.  Time  averaging  is  denoted  by  angle  brackets. 

The  variance  of  a  quantity  x,  v(x),  is  calculated  with 

v(x)  =  <x^>  -  <x>^  .  (14) 

Variances  were  calculated  for  the  horizontal  wind  speed  (u) ,  the  vertical  wind 
speed  (w) ,  fast  temperature  (T^) ,  and  one  Ly-a  sensor  (Q^) . 

The  skewness  of  temperature  and  humility  was  calculated  by 
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(15) 


S(x) 


3  3 

<x  >  -  <x> 


,3/2  • 


lv(x)] 

Structure  parameters  were  calculated  as  follows: 


2  22  2/3 

CQr  *  !‘<W  ’  -  'VV  >/r  ' 

CQT  ’  -  0A(t  «)12>/«u>t)2/3  , 


C^r  =  <(AT)2>/r2/3  , 


cl  =  <[T.(t  )  -  T.(t  +t)]2>/(<u>t)2/3 
IT  A  O  A  o 


CTQr  ‘  <<M)  ‘V<)B)>/r2/3  • 


(16) 

(17) 

(18) 

(19) 

(20) 


CTQt  “  <^TA(to)  ”  TA(to+T)1  X 


[QA(to)  '  QA(to+T)3>/(<U>T) 
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(21) 


Structure  parameters  calculated  from  sensors  separated  by  a  distance  r  are  denoted 

by  the  subscript  r;  those  calculated  with  time  delay  are  marked  with  subscript  T. 

In  formula  (16)  the  mean  difference  between  the  two  Ly-a  hygrometers  has  to  be 

included,  because  generally  <Q  >  ^  <Q  >.  The  temperature  difference  signal  AT  is 

filtered  by  a  1/60-Hz  high-pass  RC  filter,  which  guarantees  that  <T  -T  >  =  0. 

A  D 

The  vertical  fluxes  of  heat  (Qq)  and  water  vapor  (Eq)  are  calculated  by 


and 


Q  =  <wT  >  -  <w><T  > 
o  A  A 


Eo  =  <wQa>  -  <w><Qa> 


(22) 

(23) 


The  average  of  the  vertical  wind,  <rw>,  is  included  since  it  is  not  necessarily 
zero.  Apart  from  an  instrumental  offset,  <w>  will  differ  from  zero  because  of 
the  expansion  of  air  while  heated  (or  compression  if  cooled).  Moreover,  the 
sonic  anemometer  has  been  aligned  with  the  vertical,  and  is  thus  not  perpendicular 
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to  the  slightly  sloping  (1  degree)  terrain.  This  may  result  in  an  offset  in  the 
vertical  wind,  because  the  mean  wind  follows  the  terrain  and  may  have  a  vertical 
component. 


The  correlation  coefficient  r^^  of  temperature  and  humidity  was  calculated  by 


<taV  -  <TA><(V 


TQ 


[vCrA)  v(qa)] 


1/2 


(24) 


2 

The  0  of  the  10.7  and  of  the  337-pm  scintillometers  was  calculated  according 
to  equation  (13).  The  output  signal  of  the  0.94-pm  scintillometer  had  only  to  be 
averaged. 


5.  RESULTS  OF  THE  IN  SITU  MEASUREMENTS 

5.1  The  Temperature  Structure  Parameter  and  Spectrum 

2 

In  Section  4  two  methods  for  measuring  are  mentioned:  the  time-lagging 

method,  using  one  sensor,  and  the  method  with  two  spaced  sensors.  In  Fig.  8  the 

2  2 

results  of  the  two  methods  are  compared  by  plotting  the  quantity  CTl/cTr  as  a 

function  of  <u>T ,  which  Is  the  product  of  the  mean  windspeed  and  the  time  delay. 

This  product  is  equivalent  to  a  spatial  spacing.  The  fixed  distance  between  the 

2 

two  sensors  for  measuring  was  0.31  m.  In  both  examples  the  data  points 

2  2  2  2 
scatter  around  the  line  cTx/cTr  =  1  up  to  <u>i  'v  1  m,  where  CTt  drops  below  C,^. 

Probably,  at  <u>T  >  1  m  outer-scale  effects  enter,  which  means  that  the  contributing 

eddies  are  larger  than  inertial  subrange  eddies.  In  this  report,  the  structure 

-2/3 

function  multiplied  by  spacing  to  the  -2/3  power,  i.e.,  D(<u>t) [<u>t ]  is 
loosely  referred  to  as  the  structure  parameter;  strictly  speaking,  this  quantity 
is  equal  to  the  structure  parameter  only  when  the  spacing,  e.g.,  <u>x,  is  in  the 
inertial  subrange  of  spatial  sizes. 

-2 

The  measurements  at  small  <u>x  (the  lowest  value  observed  is  'VLO  m)  may 

be  contaminated  by  noise.  It  was  observed  that  the  noise  of  the  tempera ture- 

2 

measuring  equipment  was  equivalent  to  a  (0.03°C)  variance  of  the  time-lagged 

2  -2/3  2 

temperature  difference,  resulting  in  a  noise  value  of  0.2°C  m  for  at 
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2 

<u>T  *  10  '  m.  The  observed  values  of  C  at  that  spacing  were  of  the  order  of 
2  -2/3 

0.1°C  m  ,  so  a  noise  contribution  of  20%  is  likely.  An  increase  of  that 

magnitude  may  indeed  be  recognized  in  the  January  example  in  Fig.  8.  In  Fig.  9 

the  theoretical  dependence  of  the  structure  parameter  on  r/£  is  shown  in  the 

2/3  ^ 

case  of  the  Hill  spectrum.  A  decrease  of  D  (r) /r  '  is  anticipated  at  r/£  <  2. 

-2  f  ° 
Putting  £  =  0.6  x  10  m  (a  value  appropriate  to  the  left-hand  data  points  in 

0  2/3 

Fig.  8),  one  arrives  at  the  conclusion  that  D„(r)/r  should  be  nearly  constant 

-2  T 

down  to  r  ^  1.2  x  10  m.  The  observations  seem  to  be  in  accordance  with  this 
conclusion. 

2/3 

The  proportionality  of  D^(r)  with  r  is  a  direct  consequence  of  the  shape 
of  the  temperature  spectrum  in  the  inertial  subrange.  The  structure  function 
DT(r)  and  the  one-dimensional  spectrum  Ft(k)  are  related  by  the  expressions 

OO 

DT(r)  =  8tt  /  (1  -  S-}T<r)  4>t(k)k2  dK  (25) 

and  0 

OO 

F  (<.)  =471  /  <j>  (K)  K  dK  ,  (26) 

1  X  k c  1 

1 

where  ^(k)  is  the  spatial  spectrum.  These  relations  are  valid  only  if  the  tur¬ 
bulence  is  isotropic.  In  the  inertial-convective  range  (Corrsin,  1951) 

4>t(k)  =  (4tt  k2)"1  0  X  e"1/3  <~5/3  ,  (27) 

where  x  is  the  rate  of  dissipation  of  temperature  variance,  £  the  rate  of  dissipa¬ 
tion  of  turbulent  kinetic  energy,  and  0  a  constant  (8  =  0.72  according  to  Hill, 
1978b).  Combining  Eqs.  (25 )  —  (27 )  gives 

DT(r)  =  (9/10)  T(l/3)  8  XT  e"1/3  r2/3  ,  (28) 

2/3 

which  demonstrates  that  DT(r)/r  is  independent  of  r.  Let  us  now  suppose  that 
4>t(k)  ^  K^  2.  It  is  easily  found  that 
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<u>xlm) 


Figure  8. — Two  examples  of  time-lagged  observations  of  the  temperature  structure 

parameter  (c£  ),  divided  by  spaced-sensor  observations  of  the  same  quantity 
o  IT 

(C  )  versus  the  time-lagged  spacing  <u>T.  The  spacing  r  was  0.31  m. 


e  9. — The  temperature  structure  function  D^Cr)  divided  by  r  ,  and  normalized  ti 
parameter  C  versus  r/i  ,  where  i  is  the  inner  scale.  From  R.J.  Hill  (unpubl 


CO 


(29) 


DT(r)  -v  r-(l>+1)  /  (1  -  SiU)  yp  dy  . 

o  * 

The  integral  converges  for  -3<p<-l.  Then,  D^(r)  ^  r  .  The  one-dimensional 

spectrum  has  the  form 


Ft(k)  x<r  e-1/3  kp  .  (30) 

Relations  (29)  and  (30)  demonstrate  the  link  between  the  wavenumber  dependence  of 
the  spectrum,  and  the  distance  dependence  of  the  structure  function. 

The  inertial  subrange  behavior  of  the  temperature  spectrum  was  observed 

directly  by  spectral  analysis  of  the  temperature  signal.  During  five  runs,  the 

signal  was  recorded  on  an  analog  tape  recorder  and  digitized  afterwards  at  400  Hz. 

Before  digitizing,  the  signal  passed  a  low-pass  filter  with  half-power  point  at 

100  Hz  and  a  negligible  transmission  at  the  Nyquist  frequency  (200  Hz).  Each 

run  of  about  15  minutes  was  cut  into  sections  of  10  seconds  duration.  An  FFT 

procedure  was  applied  to  each  section.  Finally,  the  spectra  of  the  sections  were 

averaged  to  produce  the  spectrum  for  the  whole  run.  A  consequence  of  this  procedure 

is  that  no  information  is  obtained  for  frequencies  lower  than  0.1  Hz.  One  of  the 

temperature  spectra  is  shown  in  Fig.  10.  In  the  dimensionless  frequency  interval 

-5/3 

nz/<u>  A.  1  to  50  the  spectrum  follows  the  n  inertial  subrange  prediction.  At 
frequencies  higher  than  50  the  cutoff  of  the  electronic  filter  becomes  noticeable. 


5.2  The  Humidity  Structure  Parameter  and  Spectrum 


Results  of  the  time-lagged  and  spaced-sensors  techniques  for  measuring  the 

humidity  structure  parameter  are  shown  in  Fig.  11.  As  with  the  temperature  case, 

the  spacing  between  the  two  sensors  was  0.31  m.  For  the  February  data  in  Fig.  11, 

2  2 

it  is  seen  that  the  observed  values  of  C„  /C„  scatter  around  unity  for  <u>T 
_1  Qx  Qr 

between  10  m  and  2  m,  but  the  January  data  lie  ^15 %  above  unity.  The  difference 
between  the  two  month's  data  can  hardly  be  explained  by  instrumental  effects 
(like  the  effect  of  <Q^>  ^  <Qg>,  be  discussed  later).  Possibly  a 

anisotropic  behavior  of  the  atmosphere  is  involved,  resulting  in  larger  values 
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20  Jan  1981 
12  runs 


<u>T(m) 


2  Feb  1981 
12  runs 


2  2  2 
of  Cq  if  observed  streamwise  (Cqt)  than  observed  across-wind  (Cq^) •  The  January 

observations  in  Fig.  11  were  done  at  a  lower  wind  speed  (1-3  m  s  'S  than  those 
in  February  (3-9  ms'*'),  so  the  anisotropic  character  could  disappear  with  in¬ 
creasing  wind  speed.  Any  anisotropy  should  also  be  evident  when  comparing 
2  2 

C  /C  for  the  two  data  sets  in  Fig.  8;  however,  no  deviation  from  unity  is 

1  ir  2  2 
evident  at  inertial  range  spacings  in  Fig.  8.  Another  cause  of  might 

be  failure  of  Taylor's  hypothesis  if  the  wind  speed  is  not  constant.  The  time- 

lagged  structure  parameters  have  to  be  reduced  by  a  correction  that  depends  on 

o(u)/<u>.  The  corrections  are  the  same  as  the  ones  calculated  by  Wyngaard  and 

Clifford  (1977)  for  the  spectra.  In  the  present  case,  such  corrections  are  only 

a  few  percent  at  most. 


2  2 

For  <u>  i  >  2  m,  C  /C  drops  below  unity,  similar  to  the  behavior  of 
2  2  ^ 1  2  2 

C^/C,^.  Also  for  <u>T  <  0.1  m,  Cqr/Cqr  is  less  than  unity.  This  latter  drop-off 
at  small  spacings  is  probably  due  to  the  response  characteristics  of  the  Ly-u 
hygrometer,  which  also  affect  the  observed  humidity  spectra  (Fig.  12).  Starting 
from  dimensionless  frequency  nz/<u>  j>  0.5,  the  spectrum  follows  the  -5/3  power  law, 
up  to  nz/<u>  10;  from  then  on,  the  observed  spectral  density  is  too  low.  The 
spectral  density  is  about  half  its  expected  value  at  n  =  50  Hz  (observed  from 

_3 

five  spectra).  The  cutoff  is  as  if  the  Ly-a  had  a  time  constant  of  ^3  x  10  s. 

Using  this  time  constant  and  a  mean  wind  speed  of  4  di  s  \  an  unpublished  calculation 

2 

by  Wyngaard  shows  that  the  observed  value  of  C  is  expected  to  be  70%  of  its  true 

-2  ''r 

value  at  a  sensor  spacing  of  r  ^  6  x  10  m.  This  compares  rather  well  with  the 

-2 

observed  spacing  <u>T  ^  4  x  10  m  for  a  70%  response. 


_3 

The  "effective"  time  constant  of  3  x  10  s  is  larger  than  the  response  time 

_3 

of  the  Ly-a  as  observed  in  the  laboratory  (i.e.,  0.3  x  10  s,  see  Section  3). 

The  effect  of  space  averaging  should  be  small;  according  to  calculations  of  Andreas 
(1981)  for  a  cylinder  having  a  diameter  of  0.8  cm  and  a  length  of  1  cm  (values 
appropriate  to  the  Ly-a),  the  spectral  density  is  only  slightly  under-estimated  at 
large  nz/<u>  (see  Fig.  12).  Also,  the  effect  of  flow  distortion  was  investigated. 
For  this  purpose,  we  compared  the  spectra  of  the  temperature  and  the  humidity  of  a 
measurement  where  the  Pt-wire  was  inside  the  gap  of  the  Ly-a  hygrometer,  with  those 
of  successive  measurement  where  the  temperature  sensor  was  at  its  usual  position,  4 
cm  upstream  of  the  Ly-a.  If  there  were  an  effect  of  flow  distortion,  it  would 
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hygrometer,  with  those  of  successive  measurement  where  the  temperature  sensor 
was  at  its  usual  position,  4  cm  upstream  of  the  Ly-oi.  If  there  were  an  effect 
of  flow  distortion,  it  would  likely  be  present  in  both  the  humidity  spectrum  and 
the  temperature  spectrum.  Figure  13  shows  that  the  temperature  spectra  of  both 
situations  are  very  much  alike  (as  are  the  humidity  spectra,  Fig.  14),  which 
indicates  that  there  is  no  effect  of  flow  distortion  on  the  observed  humidity 

_3 

spectra.  Concluding,  we  cannot  offer  an  explanation  for  the  3  x  10  s  "effective" 
time  constant  of  the  Ly-ot. 


Some  remarks  have  to  be  made  on  the  stability  of  the  Ly-ci  hygrometer  relative 

2 

to  the  measurement  of  structure  parameters.  In  calculating  C  (Eq.  (16))  allowance 

2  ”r  2 

has  been  made  for  <Q  >  f  <Q  >  by  subtracting  <Q  -  Q_>  from  <(Q  -  Q  )  >.  How- 

A  d  A  d  A  o 

ever,  this  is  not  a  full  correction  because  sensor  drift  will  also  have  an  effect 

on  the  level  of  the  fluctuations.  An  analysis  (Appendix  III)  shows  that  considerable 

2 

errors  can  be  made  in  the  calculation  of  if  the  average  humidities  differ  by 

more  than  a  few  percent.  So,  it  is  necessary  to  adjust  the  calibration  constants 

of  the  Ly-a's  in  such  a  way  that  <Q  >  is  very  nearly  equal  to  <Q  >,  before  a 

series  of  runs  starts.  Once  started,  the  drift  correction  procedure  (see  Section 

3)  will  generally  maintain  <Q  >  =  <Q  >.  Still,  there  remains  a  pitfall.  The 

A  D 

drift  of  a  Ly-cx  may  exceed  the  level  of  turbulent  humidity  fluctuations  during  a 
run  (a  run  takes  about  16  minutes) .  This  happens  if  the  humidity  fluctuations 
are  very  low  (say,  less  than  0.5%  of  the  mean  humidity)  and  the  drift  rate  is 
normal  0^2 %  per  hour),  or  if  the  drift  rate  is  high  (over  50%  per  hour  was  ob¬ 
served  with  a  "dying"  detector)  and  the  humidity  fluctuations  are  moderate.  Then, 
while  the  drift  correction  procedure  is  still  capable  of  keeping  <Q  >  =  <QT1>  (if 
necessary,  the  drift  correction  procedure  can  be  adjusted  to  cope  with  a  high 
drift  rate),  it  cannot  suppress,  and  may  even  generate,  slow  fluctuations  with 

amplitude  equal  to  or  less  than  that  of  the  turbulent  fluctuations.  These  slow, 

2 

instrumental  fluctuations  may  contribute  considerably  to  C^;  in  fact,  an  in¬ 
crease  by  a  factor  of  2  and  more  has  been  noted.  This  is  a  problem  specific  to 
the  spaced-sensors  method;  the  time-lagging  method  is  insensitive  to  long-period 
fluctuations  whatever  the  source  may  be.  Hence,  the  time-lagging  technique  is 
preferred  here  over  the  spaced-sensors  technique. 


The  noise  of  the  Ly-a's  is  mostly  not  affecting  the  measuring  of  C^.  De¬ 
pending  on  the  humidity,  and  the  detector-source  separation,  the  noise  level  is 


14. — Observations  of  the  humidity  spectrum  with  the  temperature  sensor  in  the  center  of  the  Ly-a 
,ap  (dashed  curve),  and  4  cm  upwind  from  the  gap  (solid  curve).  The  straight  line  indicate  the 
5/3  inertial  subrange  power  law. 


-3  -2  -3 

in  the  range  of  10  to  10  g  m  (peak-to-peak) ,  which  results  in  a  noise  con- 

?  -7  -5  -3  2  -2/3 

tribution  to  C  of  5  x  10  to  5  x  10  (g  m  )  in  if  the  sensor  separation 

a  ^  _3  2  2/3 

is  0.31  m,  and  5  x  10  to  5  x  10  (g  ra  ;  ra  at  a  separation  of  0.01  m. 

2  _3  -32  -2/3 

Observed  values  of  are  usually  more  than  10  (g  m  )  m 


3.3  The  Temperature-Humidity  Structure  Parameter  and  Spectrum 


There  are  two  striking  aspects  of  the  plots  of  ("'j'QT/C'j'Qr  in  Fi-S*  15: 


CTQx/CTQr 


1.2  -  1.4  for  <u>r  between  0.3  and  3  m,  and  there  was  an  early  fall- 


off  of  at  small  values  of  <u>t 


Regarding  CTqt  >  C^,  only  anisotropic 


behavior  of  the  turbulence  can  be  offered  as  an  explanation;  we  cannot  explain, 

2  2 

however,  why  it  should  be  more  prominent  in  C  than  in  C^,  and  is  lacking  in  C^. 


The  drop-off  of  cTq r/^TQr  at  <U>T  <  0.2  m  was  likely  caused  by  the  joint 
effects  of  the  finite  distance  between  the  temperature  sensor  and  the  humidity 
sensor,  and  the  "effective"  response  time  of  the  humidity  sensor.  Due  to  the 
latter  effect,  an  extra  "spacing"  of  the  order  of  the  mean  wind  speed  times  the 
response  time  might  be  introduced  (for  example,  with  <u>  =  4  m  s  1  and  TLy_a  = 

3  x  10  3  s,  the  extra  spacing  is  1.2  cm;  the  instrumental  spacing  was  usually 

4  cm).  The  two  effects  would  act  in  the  same  direction  because  the  temperature 
sensor  was  positioned  upstream  (and  slightly  to  one  side)  of  the  humidity  sensor. 
Consequently,  one  would  expect  that  the  drop-off  occurs  at  higher  values  of  <u>T 
the  higher  the  wind  speed  is.  The  wind  speed  was  3  to  9  m  s  3  for  the  February 
data  in  Fig.  15,  whereas  it  was  1  to  3  m  s  3  for  the  January  data,  but  no  clear 
difference  in  the  drop-offs  of  the  two  data  sets  can  be  seen  (taking  into  account 
the  slightly  larger  values  of  C’tqj/C'tqi-  of  the  January  data  at  <u>T  ^  1  m)  .  Also, 
the  February  data,  if  split  into  a  high  wind  speed  part  and  a  low  wind  speed 
part,  does  not  show  the  expected  behavior.  Some  effect  of  the  sensor  spacing  may 
be  recognized  in  Fig.  16. 


The  temperature-humidity  cospectrum  was  calculated  for  several  runs.  An 

example  is  shown  in  Fig.  17.  We  have  to  make  a  reservation  here,  as  the  cospectral 

density  is  much  less  (by  roughly  a  factor  of  10)  than  the  cospectral  density 

implied  by  the  structure  parameter.  So,  any  conclusion  drawn  from  our  cospectra 

-5/3 

should  be  considered  tentative.  From  Fig.  17  it  is  seen  that  an  n  line  might 
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<u>x(m) 


Figure  15. — Two  examples  of  time-lagged  observations  of  the  temperature-humidity 

structure  parameter  divided  by  spaced-sensor  observations  of  the  same  quantity, 
versus  the  time-lagged  spacing.  The  dashed  line  in  the  second  example  cor¬ 
responds  to  a  n“2  dependency  of  the  temperature-humidity  cospectrum. 
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temperature  and  humidity 


Figure  17. — The  normalized  temperature-humidity  cospectrum.  The  straight  line 
corresponds  to  a  n-^/3  dependency  of  the  cospectrum. 
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be  fitted  through  the  data  points  for  nz/<u>  between  1  and  10.  At  higher  values 

of  nz/<u>  the  cospectrum  falls  off  too  rapidly,  which  can  be  ascribed  to  the 

4  cm  distance  between  the  temperature  sensor  and  humidity  sensor.  At  low  values  of 

-5/3 

nz/<u>  the  cospectrum  seems  to  decrease  more  rapidly  than  n  with  increasing  n, 
which  indicates  a  rapid  decrease  of  the  spectral  correlation  coefficient. 


These  observations  of  the  cospectrum  neither  confirm  nor  deny  the  existence 

of  an  inertial  subrange.  Other  observations  are  conflicting.  Wyngaard  et  al., 

—  S  /  3 

1978,  found  an  n  behavior  from  an  analysis  of  the  AMTEX  data;  Friehe  et  al., 
—  S/  3 

1975,  also  found  n  ,  but  they  consider  their  results  tentative.  Wesely  and 

-2.2 

Hicks  (1978)  on  the  contrary,  using  a  step-filter  technique,  found  n  '  .  A 

-5/3 

more  rapid  decrease  of  the  cospectrum  than  n 


was  also  observed  by  McBean  and 
Elliott  (1981),  but  the  0.5  m  spacing  between  their  temperature  sensor  and  humidity 
sensor  is  likely  to  have  affected  the  spectrum  beyond  nz/<u>  ^  1.  Also  the  mea¬ 
surements  by  Fairall  et  al.  (1980),  seem  to  indicate  a  more  rapid  fall-off  than  n 
they  believe  that  this  is  partly  caused  by  sensor  noise. 


5/ 


Finally,  we  may  mention  that  our  observations  of  by  time-lagging  are  in 

support  of  an  inertial  subrange  for  the  cospectrum,  because  they  show  that 
2/3 


Vr) ' r 


over  almost  a  decade  of  r. 
-2 


For  instance,  using  equations  (29)  and 

(30),  one  finds  that  an  n 

that  fits  the  observations  badly  (see  Fig.  15b). 


cospectrum  results  in  a  dependency  of  on  <u>T 


5.4  Relations  Between  Structure  Parameters  and  the  Turbulent 
Fluxes  of  Heat  and  Water  Vapor 


We  will  restrict  the  discussion  to  the  case  of  an  unstable  atmosphere  (posi¬ 
tive,  upward  flux  of  sensible  heat).  The  relations  between  structure  parameters 
and  the  fluxes  of  sensible  and  latent  heat  have  been  discussed  by  Wyngaard  and 
Clifford  (1978).  From  their  expressions  (30)  and  (31)  it  is  easily  found  that 
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-1/2 


-3/4 


<CT2)3'4' 


(31) 
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and 
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where 


E  -  a"3'4  <CT2)1/4  (C02)1/2' 

O  T  l) 


'~7t 

^Z 


(32) 


a  = 


-4/3  4  ,T  ,2/3 
2  3  (k£> 


(33) 


L  =  - 


u*  T 


kg(Qo  +  0.61  TEo/p) 


L6  = 


3  a, 
uA  T 

kg  Q„ 


(34) 


(33) 


2  -1 

Here  Q  is  the  vertical  flux  of  temperature  (K  m  s  )  and  E  that  of  water  vapor 

2°  -1  °  -1  2 
(kg  ms  ) .  The  quantity  u*  is  called  friction  velocity  (ms  ) ;  -pu^  is  the 

vertical  flux  of  momentum.  The  expression  in  equations  (31)  and  (32)  between 
parentheses  is  a  stability  correction,  with  L  the  Monin-Obukhov  length  including 
the  effect  of  humidity  on  buoyancy,  and  Lg  the  same  quantity  without  the  humidity 
effect.  The  stability  correction  in  (31)  has  been  determined  empirically  from  the 
Kansas  measurements  (Wyngaard  et  al. ,  1971),  but  no  distinction  between  L  and  Lg 
was  made  then.  Equation  (32)  has  been  proposed  in  analogy  with  (31),  and  is  in¬ 
directly  supported  by  measurements  by  Fairall  et  al.,  1980. 


-3/4  2  3/4 

We  investigated  (31)  and  (32)  by  plotting  a  (CT  )  /Qq  and 
a  ^^(CT2)1^(Cp2)^^2/EQ  versus  -z/L  (Figs.  18  and  19).  The  quantities  Qq  and  Eq 
were  measured  by  correlating  the  vertical  wind  fluctuations  with  temperature  and 
humidity  fluctuations  respectively.  In  calculating  a,  the  "Kansas"  value  of  k, 

0.35,  was  used,  and  and  Cq  were  measured  with  the  spaced-sensor  technique,  or  by 
time  lagging,  as  discussed  in  the  previous  sections.  The  friction  velocity  u*  was 
found  by  using  the  relations  of  Dyer  and  Hicks  (1970),  assuming  zq  =  0.04  m  (from 
unpublished  measurements  by  Wyngaard),  and  k  =  0.4  (Wieringa,  1980).  In  the  figures, 
the  stability  correction  as  given  in  (31)  and  (32)  is  indicated;  it  was  calculated 
putting  Lq  =  L.  (Except  for  a  few  runs,  the  difference  between  Lg  and  L  was 
negligible  anyway.)  It  is  seen  that  this  correction  is  too  large;  in  fact,  the 

_2 

data  points  do  not  exhibit  a  clear  dependence  on  stability  at  all  for  -z/L  >  2  x  10 
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Figure  19. — The  free  convection  prediction  of  the  humidity  flux  divided  by  the  observed  latent  heat  flux, 
as  a  function  of  stability.  Dashed  line  is  stability  dependence  according  to  Wyngaard  and  Clifford 
(1978). 


Also,  the  observations  do  not  seem  to  approach  the  unity  line  for  large  -z/L,  which 
is  most  clear  for  the  humidity  plot.  From  these  observations,  the  following  re¬ 
lations  are  suggested: 


where 


%  ■  “i3/4«Vl)3/4  f<- 1>  ' 


Eo  ■  “q3/4<et2)1/4  <CQ2>1/2  f(-  t>  • 


ot,j,  =  0.83  a  and  =  0.74  u,  and 


f(-  {) 


for  -  ^  >  0.02 


=  decreasing  in  an  undetermined  way  for 

0  <  -  f  <  0.02  . 

The  difference  between  a  and  a^,  may  be  non-significant  considering  the  experimental 

uncertainties,  but  it  is  believed  that  a is  definitely  somewhat  smaller  than  u . 

This  agrees  with  observations  by  Fairall  et  al.  (1980).  At  -z/L  <  0.02  there  are 

too  few  data  points,  and  their  scatter  is  too  large  to  propose  a  functional 

2  3/4 

relationship  to  -z/L.  Why  our  observations  of  (CT  )J  /Qq  differ  from  the  Kansas 
observations  for  slightly  unstable  conditions  (-z/L  0.2)  is  not  clear.  It 

should  be  kept  in  mind  that  the  Kansas  correction  may  represent  these  conditions 
poorly  because  of  the  low  number  of  observations  there. 


The  equality  of  the  stability  corrections  for  temperature  and  humidity  fluxes 

2  2  1/2 

is  more  clearly  demonstrated  in  Fig.  20,  where  the  ratio  of  (C^  /C^  )  and 

Q  / E  is  plotted  against  -z/L.  No  dependence  on  stability  can  be  detected. 

°  °  22  1/2 
Finally,  Fig.  21  gives  a  plot  of  (C„  /CL  )  versus  Q  /E  .  From  these  figures 

1  ^  0  0 

2  2  1/2 

it  can  be  concluded  that  (C^,  /C^  )  is  a  very  useful  estimate  of  the  Bowen 
ratio  (which  is  directly  proportional  to  Q^/E^),  independent  of  atmospheric 
stability. 
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2  2  1/2 

Figure  21.  —  (C^  /Cq  )  versus  Q0/E0,  which  ratio  is  directly  proportional  to  ihe 
Bowen  ratio  (=  sensible  heat  flux  divided  by  latent  heat  flux). 

6.  CONCLUSIONS 

The  operation  of  the  C07  and  the  HCN  scintillometer  was  hampered  by  insuf¬ 
ficient  stability  of  the  pyroelectric  detectors;  consequently  only  a  few  measure¬ 
ments  were  obtained  with  the  CO^  scintillometer,  and  none  with  the  HCN  system. 
Obviously,  the  application  of  pyroelectric  detectors  for  these  kinds  of  measure¬ 
ments  needs  further  research.  The  use  of  different  detectors  should  be  considered, 
e.g.,  a  HgCdTe  detector  for  the  CO^  receiver,  and  a  bolometer  for  the  HCN  re¬ 
ceiver.  These  detectors  have  to  be  cooled,  which  makes  the  construction  of  the 
receiver  units  more  complex. 
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The  LED  scintillometer  operated  satisfactorily,  although  special  attention 
had  to  be  paid  to  the  calibration  of  the  system.  It  was  found  that  the  transmitter 
aperture  was  not  completely  filled  with  light,  and  that  there  was  a  net  coherence 
across  this  aperture. 

2 

Experience  was  gained  with  Ly-a  hygrometers  as  a  tool  for  measuring  and 
C  (the  latter  quantity  involves  measuring  both  temperature  and  humidity).  Cali¬ 
bration  of  the  Ly-a  hygrometer  sometimes  drifts  too  much  to  make  such  measurements 

with  spatially  separated  sensors.  However,  by  time-lagging  observations  of 

2 

a  single  sensor,  one  can  also  make  measurements  of  (the  equivalent  for 

is  time  lagging  the  observations  of  a  temperature  and  a  humidity  sensor  that  are 

close  together),  and  this  method  is  not  sensitive  to  drift  of  the  Ly-a  hygrometer. 

Therefore,  the  time-lagging  technique  is  to  be  preferred  to  the  spatially  separated 

2 

sensors  technique  for  the  measuring  of  and  CTq. 

By  changing  the  time  lag,  the  equivalent  distance  between  the  points  of  ob¬ 
servation  (defined  as  the  product  of  the  time  lag  T  and  the  wind  speed  <u>)  may 

be  varied.  In  this  way  the  dependency  of  the  structure  parameter  on  distance 

2 

(<u>[)  was  checked.  C„  proved  to  be  independent  of  <u>T,  for  0.01  m  <  <u>T  < 

1  2 

1  m  (0.01  m  is  no  lower  limit,  but  just  the  smallest  distance  observed).  was 

constant  for  0.1  m  <  <u>x  <  2m,  and  was  constant  for  0.3  m  <  <u>T  <  3  m.  The 

drop-off  of  for  <u>T  <  0.1  m  could  not  be  explained.  The  drop-off  of 

at  <u>T  <  0. 3  m  is  likely  to  be  caused  by  the  4  cm  separation  between  the  humidity 

and  temperature  sensors.  It  is  noted  that  there  still  is  no  complete  understanding 

of  the  response  characteristics  of  a  Ly-a  hygrometer. 


The  power  spectra  of  temperature  and  humidity,  and  the  temperature-humidity 
cospectrum  were  calculated  for  several  runs  of  approximately  15  minutes  duration. 

An  inertial  subrange  was  found  for  the  power  spectra.  For  the  cospectrum,  however, 
a  definite  inertial  subrange  behavior  could  not  be  established  because  of  con¬ 
siderable  scatter  of  the  data  in  the  limited  frequency  interval  available.  From 
the  link  between  the  dependency  of  a  spectrum  on  wavenumber  on  the  one  hand,  and 
the  dependency  of  the  structure  function  on  distance  on  the  other,  it  may  be  con¬ 
cluded  that  the  cospectrum  should  have  an  inertial  subrange  as  well. 


The  relations  between  the  structure  parameters  C^,  and  ,  and  the  vertical 

fluxes  of  temperature  (Q  )  and  humidity  (E  )  for  an  unstable  atmosphere  were  in- 

o  o 

vestigated.  The  free  convection  predictions  worked  very  well,  not  only  for  very 

unstable  conditions,  but  also  down  to  -z/L  ^  0.02.  The  constant  of  proportionality 

21/2  21/2 

between  E  and  [ (C„  )  Cn  ]  is  found  to  be  about  10%  less  than  that  between 

°2  3/a  a  2  2  1/2 

Qq  and  (C^,  )  .  For  arbitrary  instability,  the  ratio  (C^  /C^  )  is  a  very 

useful  measure  of  the  Bowen  ratio  (that  is,  of  the  sensible  heat  flux  divided  by 
the  latent  heat  flux). 
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THE  C02  AND  HCN  SCINTILLOMETER  ELECTRONICS 
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APPENDIX  II 


THE  Ly-a  ELECTRONICS 
(after  Buck,  1980) 
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APPENDIX  III 


EFFECT  OF  UNEQUAL  Ly-a  HUMIDITIES  ON  THE  MEASUREMENTS  OF  C 


2 

Q 


AND  C 


TQ 


1.  The  humidity  structure  function 


The  humidity  structure  function  D^  is  calculated  with 


Dq  =  <(A-B)2>  -  <A-B>2  , 


(III-l) 


where  the  measured  humidities  are  denoted  by  A  and  B.  Now,  let  us  suppose  that 
humidity  A  is  correct,  but  B  is  off  by  a  factor  (1+6): 


B  =  (1+6)B* 

<B*>  =  <A> 

2  2 
<B*  >  =  <A  > 


(III-2) 

(III-3) 

(III-4) 


Inserting  these  expressions  in  Eq.  (III-l),  one  finds  the  observed  structure 
function  D°bS: 

D°bs  =  (1+6) [2<A2>  -  2<AB*>]  +  62[<A2>  -  <A>2]  . 

The  true  structure  function  is 


=  <(a-b*)2> 


<A-B*>2  =  2<A2>  -  2<AB*> 


so  the  ratio  of  the  observed  to  the  true  structure  parameter  is 


(III-5) 


,obs 


xtrue 


1+6  +  6' 


,true 


(III-6) 


where  v  denotes  the  variance  of  the  humidity  (=  <A2>  -  <A>2) .  The  ratio  v  /D^rue 
Q  Q  v 

generally  depends  on  the  distance  between  the  two  humidity  sensors,  and  it  can  be 

2 

quite  high — 10  or  more.  The  term  with  6  may  thus  contribute  significantly,  even 
at  a  low  value  of  6.  In  practice,  it  is  wise  to  have  6  not  larger  than  a  few  percent. 
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2. 


The  temperature- humidity  structure  parameter 


The  quantity  is  defined  by 


°TQ  '  <(VIB)<W>  ' 


(III-7) 


It  is  supposed  that  <T  -T  >  =  0  (forced  by  high  pass  filtering).  Following 

A  B 

the  same  line  of  reasoning  as  in  the  foregoing  section  we  find 


obs 

TQ 

^true 

JTQ 


<(TA-yy 

dtq 


(III-8) 


where  QD  is  the  corrected  value  of  humidity  sensor  B,  so  that  <Q  >  =  <QA>,  etc. 

B  ^ 

We  have  no  observed  values  of  <(T.-TIJ)Q_>  available.  An  upper  estimate  of  (III-8) 

ABB 

jc  9s 

is  obtained  if  we  suppose  <T.Q  >  «  <T„Q  >;  then 

A  B  B  B 


,obs 

TQ 

true 

TQ 


1  +  6 


■k 

<T  Q  > 
BVB 


=  1  +  6 


TQ 


rTQ°TgQ 

°TQ 


(III-9) 


Here  and  are  the  standard  deviations  of  the  temperature  and  humidity,  re¬ 
spectively.  The  quantities  on  the  right  in  (III-9)  all  have  been  measured.  It 
was  found  that  the  multiplier  of  6  may  be  up  to  10  with  a  spacing  of  0.31  m 
between  the  sensor  pairs  A  and  B.  Although  Eq.  (III-9)  is  surely  an  overestimate  , 
it  is  clear  that  in  the  case  of  the  measurement  of  it  is  wise  to  have  6  not 
larger  than  a  few  times  a  hundredth. 
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